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Abstract 

We present models of decaying particles for increasing the effective number of neutrinos N v after 
big bang nucleosynthesis but before the structure formation begins. We point out that our scenario 
not only solves the discrepancy between the constraints on N v from these two epochs, but also 
provides a possible answer to deeper inconsistency in the estimation of the matter power spectrum 
amplitude at small scales, represented by as, between the WMAP and some small scale matter 
power measurements such as the Lyman-a forest and weak lensing. We consider (a) saxion decay 
into two axions; (b) gravitino decay into axino and axion; (c) Dirac right-handed sneutrino decay 
into gravitino and right-handed neutrino. 
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I. INTRODUCTION 



Observations of the cosmic microwave background (CMB), galaxy clustering, weak grav- 
itational lensing and Lyman-a forest and so on strongly support the cosmological structures 
being formed in a universe described by the power-law ACDM model. Specifically, recent 
advancement in those observations enabled us to measure the matter power spectrum from 
the horizon scale down to about 1 Mpc in a very precise manner. 

However, as the data accumulate owing to the recent observations for instance by the 
WMAP and SDSS , possible tensions among different data sets are 

indicated. This is most easily seen in the term <7g, the normalization of the matter power 
spectrum at 8h~ l Mpc, where h is the Hubble parameter. Namely, the value of <7g derived 
from the WMAP three-year data is slightly lower than that derived from the latest analyses 
of the Lyman-a forest [9j, llO] , weak lensing ll|, [l2| and strong lensing ( "giant arc" ) statistics 
jl^ . The discrepancy surfaced when the WMAP data were updated from the first year data 
to the three-year data, with significant decrease in the best fit value of <7g (decreased from 
0.92 ±0.10 to 0.76li{];{jg) J fl- 
it is true that all these measurements which favor higher <jg than the WMAP3 value are 
likely to suffer more from systematic errors than the WMAP experiments, but when the 
ongoing efforts can succeed in decreasing the systematics, they would be more suitable for 
measuring <jg than the CMB experiments. Therefore, we can expect that we will obtain 
sufficient information to know whether the tension is solved by some systematics not yet 
accounted for or we have to invoke non-standard cosmology. We, in this paper, assume that 
the latter case is true and the present discrepancy between the WMAP3 and the observations 
at smaller scales is real. 

Then, what kind of non-standard ingredient do we need ? Actually, this has been already 
hinted at in the Lyman-a forest analysis of Ref. [l0]. In Ref. 10], extensive cosmological 
parameter estimation was conducted using the latest data set consists of the CMB, galaxy 
clustering and the Lyman-a; forest. They tested a wide range of cosmological models other 
than the flat ACDM model with the adiabatic power-law primordial power spectrum by 
placing constraints on the tensor mode, the running of the spectral index, massive neutrinos, 
the effective number of neutrinos, the dark energy equation of state, the curvature of the 
universe, cosmic strings and isocurvature modes. They have found that the observations 
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prefer these parameters to be consistent with the standard values except for one parameter: 
the effective number of neutrinos, N u . Remarkably, their 2a limit is N v = not 
allowing the standard value of N v = 3.0 at 2Aa 1 . 

This preference of a non-standard value of N v by the combined data of the WMAP3 
and Lyman-a forest is closely connected to the discrepancy of a 8 between these data sets. 
As shown in Ref. [k]], the larger value of N u enhances the best fit value of the small scale 
amplitude 2 . This enables the high <j% value inferred from the Lyman-a to reconcile with the 
WMAP3 data, which prefers the low erg when the standard N v = 3 is assumed. Although 
there is no detailed statistical analysis of combined data sets of the WMAP3 and the weak 
lensing allowing for the possibility of non-standard N„, it is reasonable to expect these 
observations too to be reconciled with the WMAP3 data by N v > 3. 

On the other hand, as is very well known, the value of N u greatly affects big bang 
nucleosynthesis (BBNY, especially the 4 He abundance, Y p . The analysis by Ref. [3], using 
Y p = 0.249 ± 0.009 [l7l, bij, has yielded N v = S.ltH (95% C.L.), in good agreement with 
the standard value while still allowing some room for non-standard values. For example, 
N v = 4, which can better fit the combined data of the WMAP3 and Lyman-a forest than 
N u = 3, is acceptable. However, more recent analyses favor N v = 3 (see Sec. [Ill for more 
detailed discussion). 

Having seen recent observational constraints on N v from the structure formation and 
nucleosynthesis, we will now consider how the value of N v should be in order to satisfy 
these constraints. Although the simplest choice would be to have N v ~ 4 before and after 
BBN, concerning the central values, it may be preferable to have N v = 3 during BBN and 
increase to N v = 4 - 5 well before the structure formation begins. We here note that BBN 
measures N u around the temperature T = O(MeV) while the structure formation data tell 
us N v in a more recent universe, T < 100 eV, at which the structure formation of the smallest 
observable scale (about 1 Mpc) begins 3 . In terms of the cosmological time, they respectively 
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Ref. [3] recently reexamined this issue using almost the same data set and found N u = at 95% 

C.L. Although the significance is lower than the one in Ref. lflj, tension with the standard value remains. 
Also, there is an independent analysis by Ref. [16] with a similar data set including Lyman-a which gives 
N„ — 5 ± 1, the 2a preference for JV„ > 3. 

2 In detail; they report this result in the amplitude at a smaller scale than 8ft. -1 Mpc, but similar correlation 
is expected between N w and a$- 

3 The present CMB data probe scales larger than O(10) Mpc or equivalently T < 10 eV. Meanwhile the 
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measure N v around 1 sec and after 10 s sec. Thus, the constraints on N v from BBN and the 
structure formation (the CMB, the Lyman-a forest etc.) do not necessarily coincide at face 
value in general. 

In this paper, to realize the latter possibility of increasing N u , we investigate models of 
particles which decay into radiation between BBN and the structure formation. Candidates 
would have a somewhat long lifetime of 1 to 10 s sec after which they decay "silently" , without 
destroying the light elements, into very light particles as copious as photons or neutrinos. We 
show such particles are found in supersymmetric extensions of theories proposed to solve the 
strong CP problem. Namely, we consider the following possibilities: (a) saxion decay into 
axions and (b) gravitino decay into axino and axion. We also show a candidate present in 
models with the right-handed neutrino which are attractive for explaining neutrino masses. 
In this case, we consider (c) Dirac right-handed sneutrino decay into gravitino and right- 
handed neutrino. In the next section we will give a review on the present observational status 
on as and N v and their possible tensions between different experiments. In Sec. II III we give 
details of the models and parameter space where N u is successfully increased while meeting 
cosmological constraints. Then, Sec. [TV] is devoted to our conclusions and discussion. 



II. OBSERVATIONAL TENSIONS IN a s AND N v 

In this section we give a brief review on the several different observations and analyses 
of <Tg and N v and their implications. The recent analyses of Lyman-a forest combined with 
the WMAP three-year data by two independent groups are in Refs. js] and [h]]. The earlier 
studies by the same groups with the WMAP first year data are in Refs. {20] and j^l ]. Those 
analyses used basically the same Lyman-a forest data sets. Their results seem to consistently 
show that Ug derived from the Lyman-a forest data prefers the higher value of the WMAP 
first year result rather than the three-year value, although two groups conclude that there are 
no statistically compelling evidence for inconsistency between WMAP3 and Lyman-a data. 
In Ref. 10], around two sigma discrepancy in the power-law ACDM model was reported but 



they concluded that the difference would be explained by a statistical fluctuation or unknown 
systematic errors. The analysis of Ref. |o| found weaker significance for the discrepancy and 



CMB alone does not practically constrain N v (the WMAP three-year alone limit is N v < 42 at 95% 



C.L. 



lli)- 
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concluded that the Lyman-a forest data can be in reasonable agreement with WMAP3. 
However, it is apparent from the figures of likelihood contours in Refs. (9), [h]] that the 
measurements of the small scale fluctuation amplitude by the WMAP3 alone and the Lyman- 
a forest alone are not fully consistent. 

Accumulating data for weak lensing, another efficient probe of erg, shows a similar ten- 
dency. It was first noted in the WMAP3 paper of Ref. [1] that the ground-based weak 
ensing measurement by the wide synoptic survey of the Canada-France-Hawaii Telescope 



11] favors higher values of erg ~ 0.8 - 1.0. In Ref. |l|], the likelihood contours on the f2 m -cr 8 



plane are drawn for the WMAP3 alone and the weak lensing alone but their overlapping 
region is small showing some degree of inconsistency. Moreover, a higher erg value is also 
preferred by a very recently released result of the space-based measurement by the COSMOS 
survey of the Hubble Space Telescope 12] . The agreement between the largest surveys from 



ground and space is remarkable and adds to the reliability of the weak lensing result of the 
high <j 8 . 

Regarding the fact that the numerous non-standard parameters other than N v studied in 
Ref. jjLOj] cannot solve the discrepancy, we consider that the universe with N u > 3 is a strong 
candidate for explaining both the WMAP3 and the observations which indicate high <jg. 

Now let us look at the constraints on N u from BBN. The value of N u greatly affects BBN, 
and in particular, the 4 He abundance Y p is quite sensitive to it. The analysis by Ref. flo1 |. 
using Y p = 0.249 ± 0.009 [l7, Q, has yielded N v = 3.1±H (95% C.L.), in good agreement 
with the standard value while still allowing some room for non-standard values. However, 
there are more recent analyses of Y p by several groups who give more stringent error bars: 
Y p = 0.250 ±0.004 [22], Y p = 0.2427 ± 0.0028 [23^ and Y p = 0.2516 ± 0.0011 [24]. We derive 
the constraints on N v from them using the fitting formula in Ref. [25| and the observed 
deuterium abundance D/H = (2.82 ± 0.27) x 10" 5 Q on the T]-N v plane. The 95% C.L. 
limits are respectively N u = 3.20^68, N„ = S.Oltoit and N v = 3.32+^4 (their own analysis 
in Ref. 24j] has yielded N v = 3.28 ± 0.16 (2a), using 7 Li data in addition; this is consistent 
with our calculation). Although it is beyond the scope of this paper to discuss whether their 
error bars are underestimated or not, we may conclude that three recent analyses of 4 He do 
not favor N v > 4. 

Thus we are led to consider a cosmological scenario that the effective number of neutrinos 
N u is increased from the standard value to N v > 4 during the time between BBN and the 
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structure formation 4 . In the rest of the paper, we will focus on several possible scenarios 
based on particle physics and discuss each model in detail. 



III. MODELS 

In this section we provide several models in which the effective number of neutrinos is 
increased from the standard value, N v = 3, after BBN but before the structure formation 
starts. To this end, we introduce a long-lived particle X with a lifetime tx in the range of 
Tx = O(l-10 8 ) sec. The lower bound on tx comes from the requirement that the additional 
radiation energy from the decay of X should not change the expansion rate before the 
neutrino freeze-out. This is because we do not want to change the standard BBN results, 
especially the 4 He abundance. The upper bound corresponds to the cosmic time when the 
comoving scale of about 1 Mpc enters the horizon. 

Before going to the discussion of each model, it will be useful to express the increase 
of N u in terms of the abundance and the lifetime of X. Let us assume that the decay of 
X produces very weakly interacting relativistic particles collectively denoted by R, which 
carry a fraction fn of the energy originally stored in X. The R particles increase the extra 
effective number of neutrinos by AN U : 



where px and p v denote the energy densities of X and the three species of the neutrinos, 
respectively. We define Ta as the temperature of photons when the decay rate Tx becomes 
equal to 3H (H is the Hubble parameter): 



where g* ~ 3.36 counts the relativistic degrees of freedom, and Mp = 2.4 x 10 18 GeV is the 
reduced Planck mass. The lifetime tx is related to the decay rate as tx = 1/Tx- To be 
precise, speaking, the energy density of the R particles also contributes to the right-hand side 
of Eq. d2J). Nevertheless we neglect it here, because it is sub-dominant as long as AN U < 1. 

4 Strictly speaking, the Y p analysis by Ref. 24] implies N v > 3 at 2a. However, since our models that we 
will present below are not affected by the value of N v at BBN, we assume N v = 3 at BBN for simplicity. 




(1) 




(2) 
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The standard value of the neutrino abundance is 

0.26 T d , (3) 



Pu 
S 



T=T d 



where s is the entropy density, and it should be noted here that Td denotes the temperature 
of photons, not neutrinos. Substituting Eq. (jBJ) into Eq. (PQ), we obtain 



or equivalently, 



It is clear from Eq. jSP that, to increase ]V„ by order unity, X must be produced with a 
sufficiently large abundance and its lifetime should be long enough. In order not to disturb 
the light element abundances, the decay into the standard-model particles (especially into 
the hadrons) must be sub-dominant or even forbidden, and the decay products must be 
"dark", i.e., their interaction with the visible particles should be very weak. In other words, 
any massive particles that decay into very weakly interacting and relativistic particles can 
explain the increase of N v if and only if they have right abundance and lifetime given by 
Eq. ©. 

In the following, we show three cosmologically consistent scenarios in which N v increases 
between BBN and structure formation by order unity. We consider (a) saxion decay into 
two axions; (b) gravitino decay into axino and axion; (c) Dirac right-handed sneutrino decay 
into gravitino and right-handed neutrino. We investigate each model in detail below. 

A. Saxion decay into axions 

One of the most promising solutions to the strong CP problem in quantum chromody- 
namics (QCD) is the Peccei-Quinn (PQ) mechanism [27], which involves a pseudo-Nambu- 
Goldstone boson a. the axion, associated with the spontaneous PQ symmetry breaking (for 
a review, see Ref. [281]). In a supersymmetric theory, the axion forms a supermultiplet, in- 
cluding a fermionic superpartner a, the axino, and a scalar partner s, the saxion. In general, 
the saxion acquires a mass of order 7713/2 in the presence of the supersymmetry (SUSY) 



breaking 



29. 



30] . (Here 717,3/2 is the gravitino mass.) In a class of models, the saxion mainly 



decays into a pair of the axions (i.e., fn ~ 1), and these axions contribute to the extra 
effective number of neutrinos without disturbing the BBN results 5 . However, since the 
saxion may also decay into two photons, we need to examine whether the photons produced 
do not spoil the standard cosmology. 

We consider a class of models in which the PQ symmetry is spontaneously broken by 
a single PQ scalar field $, whose vacuum expectation value (VEV) sets the scale of the 
PQ symmetry breaking scale F a = ($). Here we have assumed that the VEV of $ is real 
and positive without a loss of generality. The PQ scale F a is severely constrained from 
astrophysical and cosmological considerations as 10 10 GeV < F a < ^ _1 10 12 GeV, where 9 is 
an initial misalignment angle of the axion. It is firmly bounded from below by supernova 



cooling 
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321 ]. while the upper bound comes from the axion-overclosure 



be relaxed to some extent depending on the cosmological scenarios 
Let us express $ in terms of the saxion s and the axion a as 
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34 



imit, which can 
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36) 



$ = ^exp 



(s) 



(6) 



Expanding the saxion around its VEV as s = \/2F a + s, we obtain 



V2F a 



2 2 y/2F a 

where the third term induces the saxion decay into axions. The decay rate is given by 



(7) 
(8) 



Vis 



2a] 



1 m\ 



64tt F a 2 ' 

where m s is the saxion mass. The lifetime of the saxion then is given as 

2 



(9) 



1.3 x 10 5 sec 



(10) 



100 MeVJ V!0 12 GeV, 

Since the axion superfield $ must not have a SUSY mass, the saxion is a flat direction 
and acquires only a SUSY breaking mass of the order of the gravitino mass. Therefore, in 
the early universe, the initial position of the saxion, s» = v^l^il) naturally deviates from 



5 Note that it is model dependent whether the saxion dominantly decays into the axions. For instance, in 
a model where two PQ scalar fields <3?+ and <&_, respectively charged under the PQ symmetry +1 and 
— 1, acquire VEVs as ($+) — ($-) ~ F a , the saxion decay into axions is suppressed. 



S 



that in the vacuum (s) = y/2F a . When the Hubble parameter becomes comparable to the 
saxion mass m s , the saxion starts to oscillate around the potential minimum with an initial 
amplitude, 5s,i ~ |sj — y2F a \. There is a priori no way to determine the initial displacement 
of the saxion, 5si, but it is expected to be in the range between F a and M P . 

The saxion abundance depends on the thermal history of the universe, e.g., whether the 



reheating is completed before or after the saxion starts to oscillate [37| ■ First let us assume 
that the saxion starts to oscillate after the reheating. This is the case if the reheating 
temperature T R satisfies 

T R > 2.2xl0 8 GeVf— ^— ) 1/2 , (11) 



100 MeV 

where we have used the relativistic degrees of freedom in MSSM, g* = 228.75. The saxion 
to-entropy ratio is then given by 

Ps m 2 s (S Si ) 2 /2 3T osc 



3H 2 SC M 2 



— 6 - v (lo^v) ,/2 (lo4v) 2 (|) 2 ' < 12 > 

where the subscript "osc" denotes that the variables should be evaluated when the saxion 
starts to oscillate, i.e., H ~ m s . The saxion decays into axions, increasing the effective 
number of neutrinos AN U as 

where we have substituted Eqs. ( TlOT) and (|T2|) into Eq. ([5]), and used f R ~ 1. However, as we 
will see below, it is difficult to reconcile the constraint on T R , Eq. (TIT]) , with the gravitino 
problem. 

On the other hand, if the reheating is completed after the oscillation of the saxion com- 
mences, the saxion-to-entropy ratio is given by 
p s m 2 s (ds t ) 2 /2 3T R 



3m 2 s M 2 4 



(14) 



- 2 - 2Xl °- 8Ge Kl^)(ira) 2 (f ' • (15) 
The increase in the effective number of neutrinos is expressed as 



10 MeV/ Vl0 12 GeV/ \ F a J Vl0 6 GeV 



Thus it is possible to increase N v by order unity in this scenario. 

Now let us consider the saxion decay into two photons. The decay occurs in the DFSZ 
axion model 38[], as well as in the KSVZ (or hadronic) axion model 39] if the heavy quarks 
have U(l) em charges. To be concrete, let us consider a hadronic axion model by introducing 
the coupling of $ with the heavy quarks Q and Q as 

W = k$QQ, (17) 

where k is a coupling constant 6 . We assign the PQ charges as, e.g., $(+1), Q{— 1/2), and 
Q(— 1/2). Assuming that Q and Q furnish 5 and 5 representations of the SU (5) GUT group, 
$ couples to the standard-model gauge multiplets as 

- £ int = j d 2 6 yW®W® a + h.c, (18) 

where oti = g 2 /47r are the gauge coupling constants of the standard model, and Wa are 
chiral superfields for the gauge multiplets. Thus the saxion decays into two photons with 
the rate, 

r < s - ^ - S# < 19 > 

where k = (3/5) cos 2 8\y — 0.5 (9y/ is the weak mixing angle) and we can see that the 
branching ratio of two-photon decay is B y ~ 1.7 x 10~ 7 . The injected photons may destroy 
the light elements and change the result of BBN for m s > 40MeV, while for m s < 40MeV, 



the injected photons can never dissociate 4 He nuclei 40j. To avoid changing BBN, the 



following bounds must be satisfied [42|, |43I |: 



/p.\ I 10~ 14 GeV for 10 7 sec < r s < 10 12 sec 

By(-) << , (20) 

sJ 10- 6 -10" 14 GeV for 10 4 sec < < 10 7 sec 



and the constraints from BBN are very weak for t s < 10 4 sec. If the saxion mass exceeds 
about 1 GeV, the saxion decays into gluons with the rate 

r <* - 2 » ^ 6^4 < 2i > 



6 We assume that the PQ symmetry is broken due to the VEV of $ during inflation. Then the PQ quarks 
Q and Q are not thcrmalizcd after inflation and they do not affect the timing when the saxion starts 
oscillating. 
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which is much larger than that of Eq. ([19]) . The hadronic branching ratio is Bh = a^/ir 2 ~ 
1.4 x 10~ 3 . The bound on the saxion abundance in this case is 

/ p s \ [ l(T 13 -l(r 14 GeV for 10 4 sec < r s < 10 12 sec 



s y | 1(T 9 - 10~ 13 GeV for 1 sec < r s < 10 4 sec 

Thus, if the saxion decays into gluons, the BBN constraints on p s /s become much severer. 
In particular, for m s > 1 GeV the saxion decay cannot realize AN U = 1 due to these 
constraints. Even in the case m s < 40 MeV, the energy injection is constrained from the 
CMB. If the injected photons cannot reach chemical or kinetic equilibrium due to the small 
rate of interactions with background plasma, it leads to the distortion of the CMB black 



body spectrum which is constrained from observations 41]. Hence, the constraint comes 
from the CMB in the region m s < 40 MeV, although this does not give a severe constraint. 

Let us here comment on the thermal production of the gravitino. From Eqs. (Tl3l) and 
( Fl6l) . we can see that the light saxion mass and/or relatively high reheating temperature 
are required to obtain AN U ~ 1 as long as we stick to Ssi ~ F a . Since the saxion mass is 
considered to be of the order of the gravitino mass, the gravitino mass as well must be as 
light as 0(100) MeV, and such a light gravitino is realized in gauge-mediated SUSY breaking 



models 44j . Let us assume that the gravitino is the lightest supersymmetric partic 



If the reheating temperature is too high, the gravitino may overdose the universe 



efLSP). 
43. The 



abundance of the thermally produced gravitino is [46( (see also 431 ]) 



^- i ' 9xio " I6 ( i+ 4;j(™)' (23) 

where we have omitted the logarithmic dependence on Tr, and m~ g is the gluino mass eval- 
uated at T = Tr. For 1713/2 nig, the gravitino abundance is given as 

This should be smaller than the present upper bound on the current dark matter density, 
^dm^ 2 IS 0-12 at 95% C.L. [1]. Therefore, the thermal gravitino production sets the lower 
bound on the gravitino mass for a fixed reheating temperature. Due to this bound, there 
is no allowed region in the case that the reheating is completed before the start of saxion 
oscillations. It should be noted that this constraint cannot be alleviated even in the case of 
the axino LSP 47] . This is because, although the gravitino eventually decays into the axion 
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and the axino, the decay is too late for such a light gravitino mass. In the next subsection, 
we consider a much heavier gravitino mass, focusing on the possibility that the axion and 
the axino produced from the decay of gravitino may explain AN U ~ 1. 

The axinos, in addition to the gravitinos, are also produced by thermal scattering, and 
we should check whether the axino is overproduced. Here, since we assume that the axino 
mass is of the order of the gravitino mass, we do not have to care whether the axino is the 
LSP or not. The abundance of the thermally produced axmos is calculated as Q,Q, 

F TP * 2.0 x 10- 8 l *s(T R rH0M8/a s (T R )} \ ( F a \ 2 ( T R \ 



1.1 x 10" 4 Ml0 12 GeVy \10 6 GeV, 

or equivalently, 

OTP 7,2 ,„ Yin - 2 f a.(T R ) 3 hi[0.098/a 8 (T R )] \ / m~ a \ ( F a \ ~ 2 ( T E 
i^ft_5.5XlU ^ i.ixio-4 ) iToMevJ \1V*G<N) \l&G<N, 

(26) 

where m„ denotes the mass of the axino. Therefore, if the axino mass is too large and/or the 
reheating temperature is too high, the thermally produced axino may overdose the universe. 
Using Eqs. (fIBl) and (126]) . we derive 

From this equation, the reheating temperature should be less than 10 6 GeV as long as 
AN U < 1 and 5si > F a are assumed. Note also that the axino mass in the model ffTTl) is 
smaller than m 3 / 2 unless the PQ scalar has non- minimal coupling with the SUSY breaking 
sector. So, the upper bound on T R may be relaxed if m s <C m 3 / 2 ~ m s . 

In Fig. [TJ we summarize all the constraints discussed above. Here we have taken 6si = F a 
and T R = 10 6 GeV as reference values. The thick solid black line labeled (a) shows AN U = 1. 
Note that the region above this line corresponds to AN U > 1. The dotted blue lines (b) 
denote the astrophysical and cosmological constraints on the PQ scale, and we have set 
9 = 0.1. In order to satisfy 1 sec < r s < 10 8 sec, the combination of parameters (m s ,F a ) 
must lie in the region between two thin solid red lines (c). The constraints from BBN 
and CMB provide an upper bound on F a for a fixed m s , as represented by the dot-dashed 
green line (d). The thermally produced gravitinos exceed the current observed dark matter 
abundance if 1713/2) is smaller than the value indicated by the vertical long-dashed 

yellow line (e). For m„ = 0.01 m s , the abundance of the thermally produced axinos exceeds 

12 




FIG. 1: Constraints on the parameter space m s and F a in the saxion decay scenario with Tr = 
10 6 GeV. We have chosen 5si = F a . The lines labeled (a)-(f) are defined as follows, (a) AJV„ = 1 
on this line, (b) Lower and upper bounds on the PQ scale with 9 ~ 0.1. (c) Upper line corresponds 
to t s ~ 10 s sec, and lower line corresponds to t s ~ 1 sec. (d) BBN bounds coming from radiative 
decay for 40 MeV < m s < 1 GeV and hadronic decay for m s > 1 GeV. For m s < 40 MeV, the bound 
comes from the CMB. (e) Lower bound on m s from gravitino thermal production, (f) Lower bound 
from axino thermal production for m„ = 0.01 m s . For m„ = m s , the constraint coincides with the 
line (a) accidentally. 

the current observed dark matter abundance below the long-dashed purple line (f). For 
rria = m s , the constraint from the thermally produced axinos coincides with the line (a) 
accidentally, so it is not explicitly drawn. The region below (a) is excluded if ro„ ~ m s . We 
have found regions for ~ 1 consistent with all the constraints. For 5s i = F a , they are 
1 MeV< m s < 1 GeV, F a ~ 10 12 GeV and 10 5 GeV< T R < 10 6 GeV. Since the PQ scale F a 
is close to the upper bound coming from the axion-oveclosure limit, the axion can also play 
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F a (GeV) T R = lOMeV bs t = 10 17 GeV 




0.001 0.01 0.1 1 10 

M s (GeV) 



FIG. 2: Same as Fig. Q] except for T R = lOMeV and Ssi = 10 17 GeV. 

a role of dark matter of the universe. Moreover, dark matter may be also explained by the 
thermally produced axinos (see Eq. fl27|) ). 

On the other hand, 5si may be as large as the Planck scale. We also show the result 
for Ssi = 10 17 GeV and Tr = 10 MeV in Fig. [2} while keeping the other parameters the 
same as in Fig. [TJ Note that for such a low reheating temperature, the axion is diluted 



and the upper bound on the PQ scale F a is relaxed [36]. Since 5s{ is independent of F a 
in this case, the energy density of the saxion is also independent of F a . Hence, the BBN 
and CMB constraints is given in terms of t s and the dependence on F a is only through 
t s , as is seen from Fig. [2j For such a low reheating temperature, the axino does not give 
any meaningful constraint. For 5si ~ 10 16 -10 18 GeV, there are allowed parameter regions, 
1 MeV < m s < 1 GeV, 10 10 GeV < F a < 10 14 GeV, and a few MeV < T R < 100 MeV 7 . 



7 For the precise value of the lower bound on Tr, see Refs. 19, 50, 51, 52 |. 
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B. Gravitino decay into axino and axion 

Next we consider the gravitino decay into the axion and the axino at late times [47| . The 
axino mass is model dependent, and it can be (much) smaller than the gravitino mass [30] . 
Here, from a phenomenological point of view, we treat the axino mass as a free parameter, 
but it should be kept in mind that one may need to contrive a model that realizes a specific 
value of the axino mass, especially if it is much smaller than the gravitino mass. If the 
gravitino is the next-to-lightest supersymmetric particle (NLSP) and the axino is the LSP, 
the gravitino decays into the axion and the axino. Both the axion and the axino produced 
by the gravitino decay contribute to the extra effective number of neutrinos, so fn = 1. 

The lifetime of the gravitino is 

/ 1 VTLn In \ 

t(G -> a + a) - 



1 m l/2 
192tt Ml 



8.7 x 10 7 sec ( m3 . /2 __ ^ 



(28) 



,300 GeV/ ' 

Therefore 1713/2 must be larger than about 300 GeV in order to satisfy the requirement 
t < 10 8 sec. The needed gravitino abundance is given by Eq. (jSJ) as 

~ 8.8 x 10- 9 GeV AN„ (-^-) ^ , (29) 
s V300GeV/ ' V ' 

which is given in terms of Y as 

V 3/2 ~ 2.9 x 10~ n AN U ( m y V /2 . (30) 
3/ V 300 GeV/ V ; 

The gravitino may be produced both thermally and non-thermally. First we assume 
that gravitinos are dominantly produced by particle scatterings in thermal plasma. From 
Eq. ( J23|) . in order to obtain the gravitino abundance Eq. ( l30i) . the reheating temperature 
Tr must be as high as O(10 10 ) GeV with m- g ~ 0(1) TeV. For such a high reheating temper- 
ature, however, axinos are also efficiently produced by thermal scatterings. Their thermal 
abundance is given by Eqs. (125]) or fl26|) . Thus, for the axino abundance not to exceed the 
current observed dark matter abundance, the axino mass must be smaller than 0(1) keV. 

With such a light mass, however, its free streaming may erase the cosmological structure 
and conflict with the observation. The maximal abundance consistent with the observa- 
tional data including Lyman-a forest can be inferred from the upper bound on the HDM 



component, or the neutrino masses. According to Ref. 
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lOj, the 95% C.L. limit obtained 



from the data set including the Lyman-a forest is Yl m f < 0.17 eV which can be converted 
to Q u h 2 < 1.8 x 10~ 3 . Therefore, it is reasonable to expect that the contribution to the 
energy density of the universe from such light axino must be less than 1% of the dark mat- 
ter, in order to be consistent with the observed Lyman-a forest. This further constrains the 
axino mass down to be smaller than O(10) eV. Note that, for the axino mass lighter than 
O(10) eV, the axino abundance produced from the gravitino decay is negligibly small. 

So far we have assumed that the gravitino with the abundance Eq. (130]) is thermally 
produced. This requires a quite high reheating temperature, which limits the axino mass to 
being much smaller than the gravitino mass. Since the axino mass is generically of the order 
of the gravitino mass, such a hierarchy may pose a difficulty to build a viable axion model 
that realizes the axino mass. If the gravitino is non-thermally produced from, e.g., inflaton 



decay 



53 



54 



55 



56] (or modulus decay [57|, |58|). the tension can be greatly relaxed. The 



gravitino abundance is then dependent on the inflaton mass and VEV (0). For instance, 
in a high scale inflation model, the inflaton decays into the SUSY breaking sector, producing 
the gravitino as 



1 / / i\ \ 2 



where the precise abundance depends on the details of the SUSY breaking sector. The values 



adopted for and (0) in Eq. ( l3Ti) can be realized in e.g. a hybrid inflation model [59]. 
For such a low reheating temperature, the thermal production of the axino does not set any 
severe bound on the axino mass. In particular, note that Eq. (|25|) is not applicable for the 
reheating temperature smaller than the weak scale. Instead, the axino produced from the 
gravitino decay puts an upper bound as m a < O(100) MeV. This can be derived as follows. 
The axino abundance from the gravitino decay is 

Requiring the axino abundance to be smaller than 1% of the dark matter, we obtain an 
upper bound on the axino mass as m„ < O(100) MeV. Although the axino mass cannot be 
as large as the gravitino mass, the required hierarchy of the two is rather mild, compared 



Note that Tr ~ 10 3 GeV is naturally realized from the spontaneous decay of the inflaton through the top 
Yukawa coupling Q for to = 10 12 GeV and (<j>) = 10 15 GeV. 
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to the previous case. In a similar fashion, we can show that the LSPs produced by this 
non-thermal process cannot be the dominant component of the dark matter in the models 
described below. 

In the present model, both the axion and axino are produced from the gravitino decay 
as relativistic particles. In contrast to the axion, the axino becomes non-relativistic at some 
time depending on its mass. But it is typically well after the matter-radiation equality 
epoch, and the axino abundance amounts to only a small fraction of the energy density of 
the universe. Thus both the axion and the axino contribute to the effective number of the 
neutrinos, that is, /r = 1. The same argument is applied to the following model as well. 

So far we have neglected the saxion abundance. Since the saxion mass is roughly equal 
to the gravitino mass, the saxion decays much earlier than BBN begins (see Eq. ffTUl) ). In 
addition, if 5si is of 0(F a ), since the saxion does not dominate the universe, our arguments 
above remain unchanged. 

The final comment is that one cannot exchange the roles of the gravitino and the axino 
in the above scenario. Similar arguments show that the gravitino must be much lighter than 
the axino. This is because, as long as we require AN V ~ 1, the gravitinos produced from the 
axino are so abundant that the small scale fluctuations (> a few Mpc) would be smoothed 
out unless the gravitino mass is small enough. However, since the axion multiplet cannot 
have a SUSY mass, the axino mass cannot be much larger than the gravitino mass and the 
scenario does not seem to work. 

C. Dirac right-handed sneutrino decay into gravitino and right-handed neutrino 

The neutrino oscillation experiments have revealed that the neutrinos have finite but 
small masses. To explain the tiny neutrino masses one introduces right-handed neutrinos 
into the standard model. The right-handed neutrinos may be allowed to have large Majorana 
masses as large as GUT scale, because they are singlets with respect to the standard-model 
gauge group. However, the Majorana mass term can be forbidden by some symmetry such 
as the lepton-number symmetry. Thus, if this is the case, the neutrino mass is given by the 
Dirac mass term, and the mass of the right-handed neutrino is very light. The smallness 
of the neutrino mass is explained by the small Yukawa coupling of ~ m v j (H u ) < O(10~ 13 ), 
where m u is a neutrino mass and (H u ) denotes the VEV of the up-type Higgs. On the other 
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hand, the right-handed sneutrino acquires a mass from SUSY breaking effects. Since the 
Yukawa coupling is rather small, the lifetime of right-handed sneutrinos is very long, and 
their decay into the right-handed neutrino and the gravitino can increase N v . 
First, the lifetime of right-handed sneutrinos is given as 

-i 



Ti 



1 mf. 



" F < 1 18- mr! Al- 



3/2 iW P / (33) 



where the right-handed sneutrino mass mo R should be less than ~ 1 GeV as discussed later. 
From Eq. (J5j), the abundance of the right-handed sneutrinos should be 

* 6.8 x ID" 9 GeV AJV„ ( T {^-T . (34) 
s V500keV/ VI GeV/ 

Such large abundance of the sneutrino is unlikely to be produced by thermal scatterin gs o r 
decays of other superparticles due to the smallness of and the Yukawa coupling [60 ]. 
But, the sufficient energy density of right-handed sneutrino can be non-thermally produced 
in the form of coherent oscillations. The right-handed sneutrino can develop a large field 
value during inflation, and after the inflation ends it begins to oscillate coherently, which 



can induce a large abundance of the right-handed sneutrino [6l|, |62| . Its energy density-to 



entropy ratio Pu R /s is fixed at the end of the reheating process, 



Pur "iLfel 2 ^ 



R 



4# 2 SC M| 



4.3xHT 8 GeV( ) ( 

10 14 GeV/ V 100 GeV 



(35) 



where Um denotes the initial amplitude of the right-handed sneutrino. Here we used the 
Hubble parameter H osc at the start of the oscillations is equal to m^ R . 

As discussed above, the abundance of gravitinos produced by ur decay should be sub- 
dominant component of the dark matter. The abundance is given by 

fi3/ 2 /^9.3xl0- 4 (^) 3/2 . (36) 

It is interesting that the abundance is independent of rnz/2- in order not to significantly affect 
the observed small scale structure (> a few Mpc), m$ R must be smaller than ~ 1 GeV. With 
the constraint r < 10 8 sec, m 3 / 2 < 500 keV is also required. Note that such a hierarchical 
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mass relation may be realized in gauge-mediated SUSY breaking models with some extended 
gauge interaction which involves right-handed neutrinos and is broken at an intermediate 
scale. 

Finally, we comment on the gravitino decay into vr and Vr. This case leads to the 
same result as in the previous section, after exchanging {vr,vr) with (a, a). Therefore, for 
AN U ~ 1, 7713/2 ^ 300 GeV and < lMeV are required. However, this hierarchical 
mass relation, m 3 / 2 m i> R is unlikely in SUSY models. Hence, this case is not expected to 
explain increasing AN U . 



IV. CONCLUSIONS AND DISCUSSION 



In this paper, we present models of decaying particles for increasing the effective number 
of neutrinos N u after BBN but before the structure formation begins. In the model (a) 
where the saxion decays into two axions, broad regions are allowed. For instance, Tr can 
take a value from 10 MeV up to 10 6 GeV, depending on the initial displacement of the saxion 
and m a /m s (see Figs. [1] and [2] for details). In particular, the saxion mass needs to lie in the 
range between 1 MeV and 1 GeV, which suggests the light gravitino. In the model (b) where 
the gravitino decays into the axino and the axion, we require m 3 / 2 > 300 GeV together with 
ma ^ 10 eV or m a < 100 MeV depending on the gravitino production processes. The former 
(latter) bound is the case with the thermal (non-thermal) production. In particular, one 
needs a hierarchy between the gravitino mass and the axino mass. In the model (c) where 
the Dirac right-handed sneutrino decays into the gravitino and the right-handed neutrino, 



771,-, 



p R < 1 GeV and m 3 / 2 < 500 keV are required. This case works only for the non-thermal 
origin of the right-handed sneutrino in the form of scalar condensates. 

Such a scenario is motivated because non-standard values of N u > 3 are preferred by the 
combined data of the CMB, galaxy clustering and the Lyman-a forest 10, ll5|, 
most of the recent anal yse s of primordial 4 He abundance favor standard N u = 3 



As is discussed in Ref. 



ol whereas 
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23]. 



lOj , the preference for N v > 3 of the Lyman-a combined data stems 



from the inconsistency in the estimation of the matter power spectrum amplitude at small 
scales, represented by as, between the WMAP and the Lyman-a forest: the latter yields 
somewhat higher value of cr$. We note that such higher a$ values are also derived by other 
probes of the small scale matter power spectrum by the weak lensing III, 112J and strong 
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lensing jl3j. Thus, we would like to stress that the models proposed here can not only solve 
the discrepancy between BBN and the structure formation (the CMB, the Lyman-a forest 
and so on) but also give a possible answer to the inconsistency between the WMAP and 
small scale matter power measurements such as the Lyman-a forest and weak lensing. 

It should be emphasized that N u is increased by the "free-streaming" relativistic particles 
like massless neutrinos in our models. Our scenario of increasing N v may recall the readers 



to the scenario of "interacting" neutrinos discussed e.g. in Refs. |63l . |64| whose prediction 
includes the increase in N v after BBN by the recoupling. Even though N v changes by the 
same amount, there is a stark contrast between the free-streaming particles and interacting 
ones as regards the effects on the structure formation. The consequence is that, although 
N v can be increased in the interacting neutrino scenario, it cannot solve the problem. This 
is explicitly verified in Ref. jig . Their Fig. 5 (a) shows that the free-streaming particles can 
better fit the Lyman-a data by increasing N u from 3 but such is not the case for interacting 
particles as shown in Fig. 5 (b). 

Finally, since the discrepancy which we have addressed in this paper is about 2a level, 
further data and studies are necessary in order to see whether inconsistency exists in the 
standard cosmological model or in the interpretation of one or more observations. Recently, 



Ref. [19|| have obtained the constraint on N v from the WMAP 



QflflQ 



and the SDSS 



luminous red galaxy power spectrum [6( to be 0.9 < N u < 8.2 (95% C.L.). This is not in 
conflict with the one derived using the Lyman-a and the earlier galaxy power spectrum as 
mentioned in the Introduction, N v = 4.6^^;| [jj]], but it does not have sufficient sensitivity to 
test the need for N v > 3. Since we cannot expect the galaxy power spectrum data to increase 
significantly in near future, improvement in the reliability of the Lyman-a forest and weak 
lensing will be needed to solve the issue. We believe the ongoing works in the communities 
to understand sources of systematic errors will accomplish this task and, together with the 
uture CMB experiments (the PLANCK sensitivity for N u is forecasted to be ~ 0.2, see e.g. 



65J), this will tell us whether the scenario of increasing N v is realized in Nature. 
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